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SUMMARY 


A computer program has been developed for analyzing unresolved Mossbauer 
hyperfine spectra resulting from the existence of several local environments in 
dilute binary iron alloys. This program is an extension of an earlier program 
which had not been adapted for analyzing such unresolved compound spectra. The 
program is written in FORTRAN IV language for the Control Data CYBER 170 series 
digital computer system with network operating system (NOS) 1.1. With the pres- 
ent dimensions , the program requires approximately 36000 octal locations of 
core storage. A typical case involving two innermost coordination shells in 
which the amplitudes and the peak positions of all three components are esti- 
mated in 25 iterations requires 30 seconds on CYBER 173. The program has been 
applied to determine the effects of various near-neighbor impurity shells on 
hyperfine fields in dilute FeAl alloys. '' 


INTRODUCTION 

The effects of impurity atoms on iron Mossbauer spectra have been under 
investigation at the Langley Research Center for quite some time. (See refs. 1 
to 5.) In these studies, various types of impurity atoms in dilute concentra- 
tions (less than 6 atomic percent) were introduced into the absorber iron lat- 
tice and the Mossbauer spectra recorded. From such spectra, the effective 
values of three critical Mossbauer parameters - hyperfine field, isomer shift, 
and outer peak widths - were determined. These measurements have shown that 
(1) the effective hyperfine field and the isomer shift in iron decrease as the 
concentration of nonmagnetic impurity atoms increases and (2) the effective 
peak width shows linear increase with increasing impurity concentration. The 
effective values of Mossbauer parameters were obtained by analyzing the experi- 
mental spectrum as a single six-peak Mossbauer Spectrum, even though it was a 
composite spectrum resulting from the super imposition of spectra from several 
iron-with-near-neighbor configurations. The analysis was performed using the 
program D3290 described in an earlier report (ref. 6). Such an approach nec- 
essarily combines the effects of near-neighbor impurity atoms at various sites 
in the iron lattice. 

It would be preferable to resolve the composite iron alloy spectrum into 
single-configuration spectra so that the effects of impurity atoms in various 
configurations could be separated. Such an analysis would allow a determina- 
tion of the most important impurity sites that affect critical Mossbauer param- 
eters. It would also permit an elucidation of the mechanisms responsible for 
affecting the Mossbauer parameters. For example, it may be possible to explain 
why impurity atoms at different distances sometimes have opposite effects on 
hyperfine fields at the iron nuclei. The existing program (D3290) could not 
analyze an unresolved hyperfine spectrum into its constituents. A computer pro- 
gram that can separate individual unresolved composite peaks into constituent 
Lorentzians has now been developed and applied to the study of FeAl alloys. 



This report describes the program (M0SS2) and provides information necessary to 
use it . 
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SYMBOLS 

amplitude of absorption peak 

coefficient of second-order term in equation of parabola 
coefficient of first-order term in equation of parabola 
constant term in equation of parabola 
impurity atomic fraction 
= (Xj)2 

magnetic field at absorber nucleus 

effective magnetic field at absorber nucleus 

magnetic field at absorber nucleus in pure iron 

magnetic field at absorber nucleus in iron alloy corresponding to 
ith impurity-atom configuration 

a constant signifying effect of impurity atoms in coordination shells 
beyond 3rd near -neighbor shell 

Lorentziari profile 

number of impurity atoms in 1st near-neighbor shell 

number of impurity atoms in 2nd near-neighbor shell 

number of impurity atoms in 3rd near-neighbor shell 

number of absorption peaks in experimental spectrum 

maximum number of atoms in a near-neighbor shell 

number of total atomic configurations involved in the spectrum 

position of peak on x-axis 

probability for ith impurity-atom configuration 
probability for q impurity atoms in near-neighbor shell 
independent variable, source velocity 
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y dependent variable, counts/velocity interval 

a hyper fine field coefficient for 1st near-neighbor impurity atom 

3 hyper fine field coefficient for 2nd near-neighbor impurity atom 

Y hyper fine field coefficient for 3rd near-neighbor impurity atom 

r width of peak at half maximum amplitude 

e measurement error in L 

Xj half of peak width at half maximum amplitude for jth peak 

Subscript: 

o nominal conditions 

Superscript: 

T transpose of matrix 

Abbreviation: 

nn near neighbor 


PROBLEM DESCRIPTION 

From reference 6, a single hyperfine Mossbauer spectrum can be expressed 
as follows: 


N 

y(x) = (ax2 + bx + c) + 

j = 1 



( 1 ) 


where N is the number of absorption peaks in the spectrum. When the spectrum 
consists of a multiple set of unresolved sextets,, equation (1) must be modified 
as follows : 


y(x) 


(ax2 


6 

+ bx + c) + 

J = 1 




( 2 ) 


I 
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In this equation, each of the six hyperfine peaks is made up of n Lorentzians 
defined as follows: 


1 



n-1 




Pi 



n-1 

y' c,- 

i=0 Cj + (Pji - x)2 


n-1 

= y; Cj 

i=0 Cj + (Pji - 


where 


(3) 


Ajj^ amplitude of ith component of jth peak of hyperfine spectrum, PiAj 

Aj amplitude of jth peak for single sextet spectrum (i.e., when iron has 

no near-neighbor impurity atoms) 

Cj (Xj )2 = (r/2)2 

Pjj^ position on x-axis of ith component of jth peak 

Pj position on x-axis of jth peak of single sextet spectrum 

Pi probability for ith impurity-atom configuration 

\j half of peak width at half maximum amplitude for jth peak 

In equations (2) and (3), n designates the total number of distinct atomic 
configurations affecting the Mossbauer spectrum of absorber iron atoms. 

The analysis problem can be simplified greatly if the parabola term is 
eliminated from equation (2). This can be done by estimating the parabola coef- 
ficients with program D3290 (ref. 6). When the estimated parabola is subtracted 
from the observed spectrum, the remainder consists of six unresolved composite 
Lorentzians. Generally, the two outermost peaks are the best resolved lines in 
the spectrum and lend themselves to easy analysis. For the purpose of resolving 
these composite peaks, initial estimates of Pjj^, and Cj are needed. 

These initial estimates are obtained as follows: 

( 1 ) Initial estimate of A : 

^ji = Pi^j 

where 

Pi = p(S-,m,s,. . .) = p(L) p(m) p(s) ... (5) 
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probability 
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for 
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probability 

for 

s 


impurity atoms in 1 nn shell 
impurity atoms in 2 nn shell 
impurity atoms in 3 nn shell 


Assuming the absorber alloy to be a random binary mixture of the two types 
of atoms, probability for any given number, say q, of impurity atoms in a 
selected coordination shell, is given by the following expression: 


p(q) 



C)^max“q 


^max* cq(1 _ o^max”^ 
Q • ( ^max ” q ^ ■ 


( 6 ) 


where is maximum allowed atomic population in a near-neighbor shell 

and C is the impurity atomic fraction. Thus for any impurity-atom con- 
figuration, pj^ can be calculated from equations (5) and (6). An initial 
estimate of Aj is obtained from the experimental spectrum by analyzing 
it as a sii^le six-peak spectrum with computer program D3290 (ref. 6). 

(2) Initial estimate of Pji: 

P^i = Pi (7) 

^ Hpe 


It is shown in reference 3 that Hj_, the hyper fine field corresponding to 
the ith impurity-atom configuration, is 

f^i = %e[n + ai!- + 3m + Ys)(1 + KC)] (8) 

where Hpg represents the pure iron field obtained from the pure iron spec- 
trum using program D3290 (ref. 6), i, m, and s represent impurity-atom 
populations of 1 nn, 2 nn, and 3 nn shells, a, 3, and y represent cor- 
responding field effect coefficients, K is a constant that signifies the 
effect of all remaining coordination shells, and C is the impurity atomic 
fraction. But equation (8) cannot be used to obtain estimates of in the 

absence of information about a, 3f Y> and K. However, it has been shown 
in reference 3 that the effective magnetic field Hgff is given by 

n-1 

S PlHl 

i=0 

Heff = ^ (9) 

n-1 

Epi 

i=0 
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Initial estimates of Hgff are obtained from the experimental spectrum 
using program D3290. By combining this information with the calculated 
values of p^, initial estimates of are made. In equation (7), Hgff 

may be used as the initial estimate for Hpg; of course, Hpg can be 
obtained from the spectrum for pure iron. An initial estimate of Pj is 
obtained from the experimental spectrum by analyzing it as a single six- 
peak spectrum with program D3290. Combining the values of Hgpp, and 

Pj yields an initial estimate of Pji. 

(3) Initial estimate of Cj: An initial estimate of Cj is obtained 
from the pure iron spectrum using the program D3290 (ref. 6). It is 
assumed that Cj is independent of the impurity-atom configuration in 
the absorber lattice. 


The hyperfine field for iron alloys containing significant impurity atoms 
(i.e., greater than 5 atomic percent) is given by equation (8). However, when 
the impurity concentration is low (less than 2 atomic percent), this equation 
can be modified as follows: 


f^i = %e [( 1 + oil + 6m) ( 1 + KC)] 


( 10 ) 


This equation implies that at low impurity concentrations, only impurity popu- 
lations in the 1 nn and 2 nn shells have significant effect on the Mbssbauer 
spectrum of the alloy. Specifically, one need consider only the following 
impurity-atom configurations: (0,0), (0,1), and (1,0). Each experimental peak 

can then be assumed to be made of three peaks whose intensities are determined 
by the probabilities for the corresponding impurity-atom configurations. Thus, 
for the purpose of analyzing unresolved spectra of low-impurity-concentration 
alloys, equation (3) can be simplified as follows: 


i-j = Z 

i=0 




Cj * (Pji - x)2 


( 11 ) 


The two outer peaks (peaks 1 and 6) of the composite experimental spectrum are 
then least-squares-fitted to the sum of three Lorentzian peaks whose relative 
amplitudes are reasonably close to the estimated values of Aj^ and whose posi- 
tions are consistent with the initial estimates of Pj^. As indicated previ- 
ously, the three Lorentzians are assumed to have the same width Cj - a justifi- 
able assumption since the line width in Mbssbauer spectra is independent of 
lattice spectrum in the first approximation. 


Since equation (11) is nonlinear in several of the parameters, a nonlinear 
method was required to estimate the optimal value of the three parameters. For 
this reason the least-squares differential correction method was applied to 
this problem. A detailed description of the method is given in appendix A. 


PROGRAM DESCRIPTION 

The computer program M0SS2 is written in FORTRAN IV language for the Con- 
trol Data CYBER 170 series digital computer system with network operating system 
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(NOS) 1.1. With the present dimensions, the program requires approximately 
36000 octal locations of core storage. A typical case in which the amplitudes 
and peak positions of all three components are estimated in 25 iterations 
requires 30 seconds on CYBER 173. 

By using initial estimates for the parameters as input to the program, a 
solution is obtained. The amplitude, the position of the peak on the x-axis, 
and the square of the width of the peak at half maximum amplitude are the param- 
eters which may optionally be determined. The program allows a maximum of three 
component peaks (nine parameters) and 201 measured dependent and independent 
values. DIMENSION statements are easily changed to incorporate more measure- 
ments. In addition, the program allows input parameters describing an adjacent 
peak and a base parabola, so that these effects may be subtracted from the mea- 
sured dependent values . These are obtained from program D3290 described in 
reference 6. 


Description of FORTRAN Variables 


The following list contains a description of the significant FORTRAN vari- 
ables appearing in the program. The dimensions for each array are beside the 
variable in parentheses. The variables which are input will not be presented 
here, since they are described in a subsequent section. 


FORTRAN variable 
A(201 ,9) 
C(9,9) 

DP(9) 

DY(201) 

SUM 

TP(9) 

YP(201) 


Description 

Matrix of partial derivatives 

Matrix of coefficients for solution of simultaneous 
equations 

Right side solution of simultaneous equations 
Residuals 

Sum of squares of residuals 

Parameter values before correction by most recent least- 
squares iteration 

Dependent values computed from mathematical model 


Flow Charts and Listings of Program 

Program M0SS2 .- M0SS2 is the main program. It reads and prints the input, 
prepares for the least-squares iterations, and calls subroutine NILS. Upon 
return from NLLS, the program determines whether further cases are to be 
attempted. The flow chart and listing of M0SS2 follow. 
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1 


PROGRAM MOSS? (INPUT»0UTPUT»TAPE5-INPUT,TAP£6-00TPUT) 

COMMON A(201»9)« 0Y(201)» DP(9>» N« M, P(12},, 0€LTAP(9)> TP(9I> 

1 IPAR(9)« IT* EPS* X(201>* Y{20D* YP(20D* PA*P8*PC 
DIMENSION TITLE<8) 

5 data MMAX* NMAX/201*9/ 

C 

DO 2 J-l»201 
X(J>-0. 

2 Y(J)>0. 

10 C 

1 READ 6* (TITLE(J)* J>1*8) 

6 FORMAT(aAlO) 

IF(E0E(5i) 999*A 
A CONTINUE 

15 READ 10* 1T*H*EPS 

10 F0RMAT(2110*E20.10) 

READ 11* IPAR 

11 F0RMAT(9I5) 

DO 30^ 1-1*12 

20 READ 20* P(l) 

20 F0RMAT<F25.1A) 

30 CONTINUE 

READ 31* PA*PB*PC 

31 F0RMAT(3E16.8) 

25 READ 32* ( X ( J ) * Y ( J ) * J-1*M) 

32 F0RMAT<2E15.8t 

C COUNT NUMBER OF PARAMETERS TO FIND 

N-0 

DO 65 J-l*9 

30 IFdPAR(J) *E0. 1) N-Nd 

65 CONTINUE 

C 

PRINT A9* (TITIEUI* J-l*8) 

A9 F0RMAT(1H1//8A10//I 
35 PRINT 50* IT*H*N*EPS 

50 F0RMAT(6H IT -*I2*AH M-»IA*AH N-*I2*6H EPS-*E20.10) 

PRINT 51* IPAR 

51 FORMAT(P DELETION FLAGS... 0-DELETE* I-SOLVEP/OIIOH 
DO 60 1-1*12 

40 PRINT 55* (I*P(D) 

55 F0RMAT(5H P(*I2*1HI*1X*E18.10) 

60 CONTINUE 

PRINT 59* PA*PB*PC 

59 F0RMAT(26H PARABOLA PARAMETERS PA-*E20.8*4H PB-*E20.8*4H PC-*E20 
45 1.8) 

PRINT 61* (X(J)*J-1*M) 

PRINT 62* < Y(J )* J-1*M) 

61 F0RMAT(1H0*2H X/(6E20.8)) 

62 F0RMAT(1H0*2H Y/(6E20.8I) 

50 C 

C SUBTRACT EFFECT OF ADJACENT PEAK FROM DATA 

DO 75 J-1*M 

ADJ- P(10)*P(12)/( P(12) ♦<X( J)-P(ll))«*2 ) 

75 YU)-Y(J) -ADJ 

55 C 

CALL NLLS(NMAX*MMAX) 

GO TO 1 
999 STOP 
END 
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Subroutine NLLS .- Subroutine NLLS controls the least-squares iterations by 
counting the number of iterations, calling for function and partial-derivative 
evaluations, solving the simultaneous equations, and correcting the parameters 
being determined. Also, this subroutine prints the results of each iteration. 
The flow chart and listing of NLLS follow. 
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1 

9 


10 


19 


20 


29 


30 


39 


40 


49 


90 


99 


60 


SUBROUTINE NLLS ( MN4X, HNAX I 
C 

C NON-LINEAR LEAST SQUARES IITERATEO) 

C 

COHNON A(201>9)» OY(201)> Dl>f9)« N> N, P(12)> DELTAP(9)> TP(9)» 

1 IPAR(9)> IT* EPS* X(201l* YI201)* YP(20D* PA*PB*PC 
DIMENSION C(9*9)* IPIV(9)* IWK(IB) 

DIMENSION ABC(9) 

DATA (ABC(J)*J>1*9)/10HA0 *10H PO *10H CO *10H 

lAl *10H PI *10H Cl *10H A2 *10H P2 

210H C2 / 

DATA IP*MMN/1*1/ 

C 

DO 100 KK-1*IT 
PRINT 11* KK 

n F0RMAT(1H1*40X*4ITERATI0N *14) 

CALL FUNC 
SUM>0. 

DO 31 J-1*M 

31 SUM > SUM * 0YUIPP2 

PRINT 32*SUM 

32 FORMATdHO** SUM-OF-SOUARES 4E20.7) 

DO 20 I*1*N 

C FORM OP > A(TRANSPOSE)«OY 

0P<I)-0. 

DO 14 L-1*M 

14 DPdI-DPm * A(L*n*OY(L) 

C FORM C- A(TRANSPOSE)*A 

DO 19 J-1*N 
Cd*J)-0. 

DO 19 L«1*M 

19 C(I*J>- Cd*J) «A(L*T)4A(L*J) 

20 CONTINUE 

CALL MATINV(NMAX*N*C*MMM*DP*IP*DET*ISC*IPIV*IWK) 

IFtDET .NE. 0.) SO TO 103 
PRINT 101* OET 

101 F0RMAT<9H DET-*E15.7** STOP ♦) 

STOP 

103 CONTINUE 

PRINT 102* (C( J*J)* J"1*N) 

102 FORMATdOH VARIANCES*6E16.7/3E16.7) 

DO 110 J-1*N 

SUM - SQRT( C(J*JI) 

DO 110 L-1*N 
C<J*U- C(J*L)7SUM 
C(L*J)> C(L*J)/SUM 

110 CONTINUE 
PRINT 111 

111 FORMATdH ** CORRELATIONS 
DO 120 J-1*N 

PRINT 112* (C(J*L)*L-1*N) 

120 CONTINUE 

112 FORMATdH *9E14.7) 

C 

C LINE-UP CORRECTIONS 
K-1 

DO 10 J>1*9 
IFdPAR(J)) 9,9*9 
9 DELTAP(J)-OP(K) 

K-K+1 
60 TO 10 
9 DELTAP(J)«0. 

10 CONTINUE 
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69 


70 


79 


60 


DO 470 1-1*9 
OP(I)-DELTAP(n 
470 P(U-TP(ri* opm 
PRINT 920 

920 F0RNAT(//«1X*4ELEHENT NANE*«14X*4aL0 VALUE4*21X**1NCRENENT4> 
.21X*4NEW VALUE**/) 

PRINT 990* (ABC(J)*TP(J)*DP(JI*P(J)*J>1*9) 

990 FORMAT (1H0*(A10*8X*E21. 14*9X*E21. 14*9X*E21. 14 ) ) 

C ■ 

C IF ALL PARAMETER CORRECTIONS LESS THAN EPS STOP 

DO 600 J-1*N 

IF( ABS(OPIJ)) .6E. EPS) 60 TO 999 

600 CONTINUE 
PRINT 601*EPS 

601 F0RNAT(27H ALL CORRECTIONS LESS THAN *E19.7*7H STOP.) 

RETURN 

999 CONTINUE 
100 CONTINUE 
RETURN 
END 


Subroutine FUNC .- Subroutine FUNC uses the parameters to evaluate the pre- 
dicted dependent variable and all the partial derivatives and evaluates the 
residuals. The matrix of partial derivatives is compressed according to which 
parameters are not to be determined. The flow chart and listing of FUNC follow. 


G5ED 


Evaluate dependent variable, 
partial derivatives, and 
residual at each 
value of X 



Delete c( 
correspond! m 
not being 

jlumns of A 
3 to parameters 
determined 
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5 


10 


15 


20 


25 


30 


35 


AO 


SUBROUTINE FUNC 

CONNON AC201»9)» DY(201)» DP(9I# N# N> P(12)> DELTAP(9)» TP(9I» 
1 IPAR(9)> IT* EPS* X(20D* Y(20D* YP(20D* PA*PB*PC 
C 

C EVALUATE FUNCTION AND PARTIALS 

C 

00 60 1-1*9 
60 TP(I)-Pm 
C 

DO 10 J-1*H 
YPU)-0. 

00 7 K-l*3 

L-3*K 

L2-L-2 

DO- P(L-1)-X(J) 

C PARTIAL WITH RESPECT TO AIK) 

A(J*L2)- P(L )/(P(L)«00«*2) 

YP(J)- YP(J)« P(L2)*A(J*L2) 

C PARTIAL WITH RESPECT TO PIK) 

A( J*L2^D- -2.P00«P(L2)*A<J*L2)/(P(L)400**2) 

C PARTIAL WITH RESPECT TO CIK) 

A(J*L) ■ P(L2)*00**2/I (PIL )«00P«2)'*«2 ) 

7 CONTINUE 

C ADO PARABOLA TO FUNCTION 

TPIJ)- YP(J)+ PA*X(JI*A2 4P6*X(J) 4PC 
C RESIDUAL \ 

OYU)- YU)-YP(J) 

10 CONTINUE 
C 

C DELETE PARAMETERS 
C 

K-0 

00 40 1-1*9 

IF( IPAR(I) ) 40*40*20 
20 CONTINUE 

C SAVE COLUMN I BY MOVING TO COLUMN K 

K-K+l 

DO 29 J-1*M 
29 AU*K)- A(J*I) 

40 CONTINUE 
RETURN 
END 


Subroutine MATIHV .- Subroutine MATINV is described in detail in appendix B, 
Subroutine NLLS uses this subroutine to invert and solve the matrix of normal 
equations . 


PROGRAM USAGE 
Input 

All the data are input to M0SS2 by using FORTRAN READ statements. For all 
the input variables that follow, the format used is shown in parentheses. 
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Card 

FORTRAN 

variable 

Description 

Format 

1 

TITLE 


Case identification 

(8A10) 

2 

IT 


Maximum number of iterations 

(2110, E20 


M 


Number of data- points 



EPS 


Stopping criteria; the program stops 





iterating on this case if the abso- 
lute value of each correction is 
less than EPS 


3 

IPAR 


Deletion flag: for 1 5 I 5 9, 

(915) 




IPAR(I) =0 Do not solve for Ith 





parameter 





= 1 Solve for Ith parameter 

4 

P(1) 


Amplitude of first component peak 

(F25. 14) 

5 

P(2) 


Position of first component peak 

(F25.14) 

6 

P(3) 


Square of half-width of first compo- 

(F25.14) 




nent peak 


7 to 12 

P(4) to 

P(9) 

Amplitude, position, and square of 

(F25.14) 




half-width for second and third com- 
ponent peaks (same order as cards 4 
to 6) 


13 to 15 

P(10) to P(12) 

Amplitude, position, and square of 

(F25. 14) 




half-width for adjacent peak (same 
order as cards 4 to 6) 


16 

PA 


Coefficients of parabola, 

(3E16.8) 


PB 


ax^ + bx + c 



PC 




Next M 

X 


Independent variable x 

(2E15.8) 

cards 

Y 


Dependent variable y 



Output 

An example of the output is shown in appendix C. Explanations of output 
headings are as follows; 

Heading Description 

Case identification Identification given in TITLE 

IT Input value of IT 
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Heading 


Description 


M 


Input value of M 




N 


. Number of parameters 

to 

be 

determined 

EPS 


Input 

value of EPS 




DELETION 

FLAGS . . . 

Input 

values in IPAR(1) 

to 

IPAR(9) 

P(1) to 

P(12) 

Input 

values in P(1) 

to 

P(12) 

PARABOLA 

PARAMETERS 






PA 


Input 

value of PA 




PB 


Input 

value of PB 




PC 


Input 

value of PC 




X..;.. 


Input 

values of X 




Y 


Input 

values of Y 





The following information is printed for each iteration: 


Heading 

ITERATION 


Description 
Iteration number 


SUM-OF -SQUARES 

VARIANCES 

CORRELATIONS 
ELEMENT NAME 
OLD VALUE ; 

INCREMENT 


NEW VALUE 


Value of the sum of squared residuals for 
the current iteration 

Value of the parameter variances given by 
the least-squares procedure 

Correlation matrix for the parameters 

Names of the individual parameters 

Parameter values at the beginning of the 
least-squares iteration 

Parameter corrections obtained by least- 
squares iteration if parameter was to 
be determined; if parameter was not to 
be determined, a zero is printed 

Parameter values after parameters are 
updated with correction values 


If the program stops iterating because all the corrections fulfill the 
stopping criteria, the following message is printed: 


ALL CORRECTIONS LESS THAN (Value of EPS) 


STOP. 



The program stops and goes on to the next case if either the stopping criteria 
is attained or the maximum number of iterations have been attempted. 


Applications ■ 

Mossbauer spectra of FeAl alloys containing 0.5, 1.0, 1.5, and 2.0 atomic 
percent aluminum were studied with a conventional constant-acceleration trans- 
mission Mossbauer spectrometer. The two outer peaks were analyzed into three 
constituent Lorentzians at each impurity concentration. From the positions of 
the three Lorentzian peaks, the hyperfine field values in iron resulting from 
1 nn and 2 nn aluminum atoms were determined. The results are summarized in 
tables I and II and are shown for the 2.0-percent aluminum alloy in figures 1 
and 2. These results are in reasonably good agreement with the reported values 
(refs. 7 to 9) . 
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TABLE I.- SUMMARY OF CONSTITUENT LORENTZIAN PARAMETERS FOR PEAKS 1 AND 6 IN FeAl 


MOSSBAUER SPECTRA AT SEVERAL ALUMINUM CONCENTRATIONS 


Aluminum 
concentration , 
atomic percent 

Amplitudes of three 
Lorentzian peaks, 
peak counts 

Positions of three Lorentzian peaks 
on x-axis, mm/sec 





Absorption peak 1 


,0 

Al,l 

^1,2 

Pl,0 

P 1 .I 

. Pi, 2 

0.5 

4723 ± 63 

00 

1+ 

142 ± 63 

-4.8167 t 0.0024 

-4.3644 ± 0.0604 

-4.5570 ± 0.0864 

1.0 

4544 ± 66 

364 ± 66 

273 ± 66 

-4.7960 ± 0.0026 

-4.4115 ± 0.0338 

- 4.5870 ± 0.0476 

1.5 

4305 ± 66 

517 ± 66 

387 ± 66 

-4.7902 + 0.0029 

-4.4209 ± 0.0249 

- 4.6812 ± 0.0358 

2.0 

4388 ± 66 

702 + 66 

526 ± 66 

-4.7942 ± 0.0029 

-4.3967 ± 0.0194 

-4.6035 ± 0.0267 


Absorption peak 6 



A6,1 

*6,2 

P6,0 

P6,l 

P6,2 

0.5 

5005 ± 55 

200 + 55 

150 + 55 

+4.2396 ± 0.0016 

+ 3.8770 ± 0.0429 

+4.1251 ± 0.0621 

1 .0 

4688 ± 48 

375 + 48 

281 ± 48 

+4.2337 ± 0.0020 

+ 3.9353 ± 0.0213 

+ 4.1702 ± 0.0369 

1.5 

4696 ± 48 

564 ± 48 

423 ± 48 

+ 4.2349 ± 0.0017 

+ 3.9001 + 0.0147 

+4.0607 ± 0.0204 

2.0 

4656 + 45 

745 + 45 

559 ± 45 

+ 4.2314 ± 0.0017 

+ 3.9058 ± 0.0109 

+4.0412 ± 0.0149 
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TABLE II.- SUMMARY OF MEASURED VALUES OF AH FOR 1 nn AND 2 nn 


ALUMINUM IMPURITY ATOMS IN FeAl ALLOYS 


Aluminum 
concentration , 
atomic percent 

Ah, percent reduction, for - 

1 nn shell 

2 nn shell 

0 

0 

0 

.5 

8.0 ± 1.3 

3.9 ± 1.8 

1 .0 

7.3 ± 0.8 

2.9 ± 2.0 

1.5 

7.6 ± 0.6 

2.9 ± 0.9 

2.0 

8.0 ± 0.4 

4.3 ± 0.6 

Average . . . 

S 7.7 ± 0.8 percent 

^ 3.5 ± 1.3 percent 


^These values should be compared with previously 
reported reductions of 8.0 and 4.0 percent, respectively 
(refs. 7 to 9). 




COUNTS/VEL(X)TY INTERVAL 


100 


Z L(o,ci+L(i.o)+L(o,i) 


3200 


2400 


SOURCE VELOCITY, mm/sec 

Figure 1 Illustrative analysis of a composite Mossbauer hyperfine peak 
into three constituent peaks in FeAl alloy with 2.0 atomic percent 
aluminum. 




SOURCE VELOCITY, mm/sec 


Figure 


- Comparison between the computed and the experimentally observed 
spectra in FeAl alloy with 2.0 atomic percent aluminum. 





APPENDIX A 


LEAST-SQUARES SOLUTION 

The solution of equation (3) contains 3n parameters, which under theoretical 
conditions are constant. These parameters are AqjPo.Co,. . • iAn_i >Pn-1 >*^n-l 
where Cq = Ci = . . . = Cn_i = r2/4. 

Associated with L in equation (3) is some measurement error e. Since 
this error exists, a computational method is needed which yields the best pos- 
sible results with all the information available. The method of least squares, 
which is described subsequently, uses a minimum error criterion and has been 
used in the analysis of these Mossbauer data. 

In general, equation (3) with the associated measurement errors can be 
written as 


Lj^ = F(x , A q , Pq » Cq > • • • > ^n-1 > Pn-1 » ^n-1 ) "*■ ^i (A1) 

where F is a nonlinear function of the parameters, and in order to find an 
estimate to the parameters of the function, it must be linearized. Expanding 
equation (3) in a Taylor series about a nominal set of parameters and dropping 
the higher order terms result in the following linear approximation: 

ALi = Li - Li^o 

= bijCAo - Aq^q) + ^i,2(P0 - Pq,o) + bi^3(Co - Cq^q) + • • . 


+ *3i,3n-2(An-i - *^i,3n-l(^n-1 “ Pn-1,o^ 

'^i,3n^*^n-1 “ *^n-1,o^ ^i 


where 




Aq 


o 





3Fi 

Co 


o 


^i,3n-2 


_ 3Fi 


^n-1 


bi,3n-1 


3Fi 

Pn-1 


o 


h 

^i,3n = — 


Cn-1 
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If ttj (j = 1, 2, . . 3n) represent the 3n parameters and 

Aa-| = Aq - Aq,o ^2 = Pq - P0,o ' ^“3 = - Cq,o 

^°^3n-2 “ ^n-1 “ An-1,o ^°‘3n-1 = Pn-1 “ Pn-1,o ^*^3n = ^n-1 ” *^n-1,o 

equation (A2) can be put into the following form: 

3n 

ALj^ = ^i, j ^i (A3) 

j=1 ’ 

For further considerations, the linear equation (A3) corresponding to the ith 
observation is expressed in matrix notation as 

= Bi Aa + e^ (A4) 

where 

Zi = Ayi Bi = •] ,bi^2> • • •>*^i,3nj 



Then for k (k ^ 3n) observations, there are k matrix equations of the form 
of equation (A4) which may be written 

z = B Aa + e (A5) 
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where 






ei 

Z2 


B2 


®2 

• 

B = 

■ 

e = 

• 

_zk_ 


3c_ 


6k_ 


The problem may be restated ; Given z and S, find the best estimate Aa 
for Aa . 

Th^ best estimate Aa is the value of Aa which minimizes the sum of 
the squares of the residuals e'^e where 

e'^i = (z - B Aa)'^(z - B Aa) (A6) 

In order to minimize equation (A6), the first variation 5 with respect to Aa 
must vanish; that is, 

6(eTi] = (s[(z - B Aa)'^(z - B Aa)] = -2(£T - AaTBT)B 6(Aa) = 0 (A?) 

Since 6(Aa) i o, equation (A7) can be satisfied if 
(z'^ - AaT bT)B = 0 


S'^B Aa = b'^z 


(A8) 


Solving for the estimate of Aa in equation (A8) gives 

^ = (bTb)“^bT£ ' (A9) 

A second necessary condition for equation (A6) to be a minimum is that the 
second variation with respect to Aa be positive definite. Upon examination, 
the second variation is 

6.2(eTi) = 2 6(Aa)T gTg 5(Aa) 

jihich is positive definite. Therefore, equation (A9) is a valid expression for 
Aa. 


^ Since equation (A9) is based on a linear approximation with nominal 
Aa can be used to find the^best estimates ctj. With the relationship ’ 

a = aQ + Aa, the value pf Aa which minimized equation (A6) leads to a new 
nominal • This process implies an iterative procedure which 
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continues until Act -»■ 0 and the value of that leads to this result is 

the best estimate Sj for aj. 


Error Analysis 


Associated with the least -squares^solution is the determination of the 
accuracy of the parameters ol. When M 0 the best estimate a* of a 

is obtained. Associated with a* is the matrix . Multiplying 


[bTb] 


- 1 * 


by the predicted variance 


where 


a2 _ i'^e 

(k - 3n) 


yields the covariance matrix V of the estimated parameters, 

V = [BTgJ"^*a2 

By examining the diagonal elements of V, estimates of the error in a can be 
obtained. That is, the square roots of the diagonal elements of V are the 
estimates of the error bounds of a. 


Partial Derivatives 


For equation (A1), 
n-1 




F = _ 

, j=0 Cj + (Pj - x)“^ 

The ^partial derivatives of F with respect to each parameter are as follows: 
9F Cj 


3A 


j Cj + (Pj - ,x)2 


3F = 

3Pi /r, 


J [Cj + (Pj - x) 
^ ~ 


(j=0, 1, . . .,n-1) 


fcj + (Pj - X)2]^ 
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LANGLEY LIBRARY SUBROUTINE MATINV 

Language ; FORTRAN 

Purpose ; MATINV solves the matrix equation AX = B, where A is a square coef- 
ficient matrix and B is a matrix of constant vectors. The solution to a 
set of simultaneous equations, the matrix inverse, and the determinant may be 
obtained . 

Use ; CALL MATINV(MAX,N, A, M, B, lOP, DETERM, ISCALE, IPIVOT, IWK) 

MAX The maximum order of A as stated in the DIMENSION statement 

of the calling program 

N The order of A; 1 1 N 1 MAX 

A A two-dimensional array of coefficients. On return to the 

calling program, A“^ is stored in A. 

M The number of column vectors in B. On return to the calling 

program, X is stored in B if M > 0; for M = 0, the sub- 
routine is used only for inversion. 

B A two-dimensional array of the constant vectors B. On return to 

the calling program, X is stored in B. 

lOP Option to compute the determinant: 

0 Compute the determinant. 

1 Do not compute the determinant . 

DETERM Gives the value of the determinant by the formula 

Det(a) = lo(IOOxISCALE) ( determ) when lOP = 0. For lOP = 1 , 
the determinant is set to 1 . For a singular matrix and 
lOP =0 or lOP = 1, the determinant is set to zero. 

ISCALE A scale factor computed by the subroutine to keep the results of 

computation within the floating-point word size of the computer 

IPIVOT A one-dimensional array of temporary storage used by the 

subroutine 

IWK A two-dimensional array of temporary storage used by the 

subroutine 

Restrictions ; Arrays A, B, IPIVOT, and INDEX have variable dimensions in the 
subroutine. The maximum size of these arrays must be specified in a DIMENSION 
statement of the calling program as A(MAX,MAX), B(MAX,M), IPIVOT(MAX) , and 
IWK(MAX,2). The original matrices A and B are destroyed. They must be saved 
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by the user if there is further need for them. The determinant is set to zero 
for a singular matrix. 

Method : Jordan's method is used to reduce a matrix A to the identity matrix I 

through a succession of elementary transformations 5-^, ^n-1« • • S-i- If 

these transformations are simultaneously applied to I and to a matrix B of con- 
stant vectors, the results are A“^ and X where AX = B. Each transformation 
is selected so that the largest element is used in the pivotal position. (See 
ref. (a).) 

Accuracy ; Total pivotal strategy is used to minimize the rounding errors; how- 
ever, the accuracy of the final results depends upon how well-conditioned the 
original matrix is. A return with DETERM i 0 does not guarantee accuracy 
in the solutions of inverse. 

Reference ; (a) Fox, L.: An Introduction to Numerical Linear Algebra. Oxford 

Univ. Press, 1965. 

Storage : 516 octal locations 

Subroutine date ; January 1, 1975 
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